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Abstract 

j^-^ ' We develop a systematic approach to deriving addition theorems for, and some other bilocal sums of, 

spin spherical harmonics. In this first part we establish some necessary technical results. We discuss the 
factorization of orbital and spin degrees of freedom in certain products of Clebsch-Gordan coefficients, 
and obtain general explicit results for the matrix elements in configuration space of tensor products 
of arbitrary rank of the position and angular-momentum operators. These results are the basis of the 
addition theorems for spin spherical harmonics obtained in part II. 



?H ■ 1 Introduction 

The importance of the representation theory of the three-dimensional rotation group [IHZ] in the study of 
all natural quantum systems hardly needs to be mentioned. That central role is due to the fact that the four 
^ . fundamental interactions of nature, and their effective interactions relevant to nuclear, atomic and molecular 

physics, are all invariant under the rotation group. It is also for that reason that partial-wave expansions 
. are essential tools in the analysis of both classical and quantum scattering processes. 

' The algebraic aspects of the computation of partial-wave expansions with spin states are best codified 

by the addition theorems for spin spherical harmonics. Spin-s spherical harmonics Y^j {r) [Hin] are the 
angular-momentum eigenfunctions relevant to the description of spin-s particles subject to spin-dependent 
central interactions As such, they are of interest in their own right beyond their above-mentioned specific 
application to partial-wave expansions that constitutes our main motivation for their study. In this paper, 
together with its second part [lOj (hereafter referred to as II), we develop a systematic framework to derive 
addition theorems for spin spherical harmonics, and employ it to obtain those addition theorems for low 
spins, < s', s < 3/2, and for arbitrary integer spin s' if s = 0. 

In this first part we establish the preliminary results that will be used as basic building blocks for 
constructing addition theorems for spin spherical harmonics in II. We consider first the factorization of 
orbital and spin degrees of freedom in products of Clebsch-Gordan (henceforth CG) coefficients of the 
form (^', t'^] s', s^Ij, jz){i, iz'yS, Sz\j, jz)- We obtain a general factorization result, make its tensor and spinor 
structure explicit, and discuss the particular cases relevant to spins 1/2, 1 and 3/2. Such products of CG 
coefficients occur in the computation of physical quantities, such as scattering amplitudes or other matrix 
elements, related to transitions among states with definite orbital and spin quantum numbers. Thus, these 
results may be of interest independently of the addition theorems considered here and in II. We study also the 
matrix elements of arbitrary tensor products of orbital operators between angular-momentum projections 
of position eigenstates, for which we obtain completely general results. These are the other main ingredient 
in our derivation of addition theorems for spin spherical harmonics in II. Those matrix elements are closely 
related to bilocal spherical harmonics and other bilocal sums of ordinary spherical harmonics, as discussed 
in detail in II, so we expect that these results and the approach used to derive them should be applicable in 
other contexts as well. 



o 
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^Familiar examples of spin-s spherical harmonics are ordinary scalar spherical harmonics, Y^^ (r) = Yu^lj), and vector 
spherical harmonics [ll|6l|8l|9] , Yjj (x)- 
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The outline of the paper is as fohows. In the next section we discuss the factorization of orbital and spin 
degrees of freedom in products of CG coefficients. In sectionOwe introduce the angular-momentum projector 
operator and obtain general expressions for its matrix elements with orbital tensor operators. In appendix 
l^we state our notation and conventions. In appendix IB] we give the definitions and main properties of the 
standard bases of irreducible tensors and spinors, of arbitrary spin, used as spin wavefunctions throughout 
the paper. A detailed list of general results for reduced matrix elements of arbitrary tensor products of 
orbital and spin operators, needed in the main sections of the paper, are given in appendix [C] Finally, 
appendix ID] gathers some ancillary calculations needed in sect. [3] 



2 Factorization of orbital and spin dependence 

We begin by considering products of CG coefficients of the form 

i 

with £ integer and s', s integer or half-integer. We will usually write S{j,£' s' , s), omitting the last 
arguments for brevity. For the purpose of obtaining addition theorems for spin spherical harmonics we need 
to rewrite S{j, £' , £, s' , s) as a sum of terms with a completely factorized dependence on orbital and spin 
angular-momentum projections, and to make explicit the tensor structure of those terms. The first goal is 
achieved by making use of relation (jA.l[) to write (see appendix |X] for our notation and conventions) , 

su,£',£,s',s) = i-ir''+-'^i-iY+''-^^^i^^= (2A+i)(^'^ i 

X (£,4;A,A4|^',4)(s,s,;A,As,|s',4) , 
A,„in = max{|A^|,|As|} , A^,^^ ^ mm{£' + £, s' + s} , A£, + As,=0. 

In order to make the tensor structure of each term on the r.h.s. explicit, we use the Wigner-Eckart (henceforth 
WE) theorem [IKSllH] to write the CG coefficients in terms of appropriate irreducible tensor operators, 



S{j,£',£,s',s) = Ci!^{£',£'^\e''---'^{A£,)?'' ...?'^^'^U\^'\+^ ...V^\£,£,) 

A=A„i„ 

X (s',4|e''i-''^(A4)*r''i ...r''i^=i5''i^=i+i ...5''^|s,s^) , 

Ci!^ = (-1)'^'^ (-1)^+-'-^-^= ^^' + ^ (2 A + 1) ( ^' ^ U 

V y V 7 y(2FTT)(27TT) J s J 



(3) 



where Amin.max are as in ([2]), and our notation for the standard basis tensors - and for irreducible 
tensor operators is explained in appendix |B] The operators appearing in ([3]) are the orbital L and spin S 
angular momenta, the position versor r = f/r, and the spin-transition operator T defined here as the spin 
analog of r (see sect. 12.41 below for a detailed discussion of this operator). The reduced matrix elements 
appearing in ([3]) are explicitly given for arbitrary values of their parameters in appendix [C] 

Whereas ([3]) is completely general, it is not yet explicit enough for the purpose of deriving addition 
theorems for spin spherical harmonics. From a practical point of view, we need only consider low values of 
s', s since quantum states with high values of spin occur infrequently in nature. In this paper we restrict 
ourselves to s' , s = 0, 1/2, 1, 3/2, and |As| < 1. For those spins the coefficients C|,/j^ can be drastically 
simplified by exploiting the fact that j — £' and j — £ can take only a small set of values. The spin matrix 
elements in ([3]) can be compactly evaluated in terms of standard tensors and spinors, and spin matrices. 
(The orbital matrix elements could, in principle, also be expressed in that way, but it would be impractical 
because £\ £ are not bounded above thus requiring tensor wave functions of arbitrarily large rank. We 
evaluate those matrix elements with a completely different technique below in section |3l) Clearly, we must 
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discuss separately the different values of s' and s, as those will determine the kind of spin spherical harmonics 
involved in the addition theorem. For fixed s' and s, different values of Ai correspond to different addition 
theorems, so we must consider those cases separately as well. Furthermore, both in addition theorems for 
spin spherical harmonics and in partial wave expansions, terms with different |A£| have different tensor and 
spin structure. 

2.1 s' = l/2 = s 

For s' = 1/2 = s there are two possible values of \A£\ = 0, 1. 

A£ = In the case s = 1/2, j can take only two values and ([3]) reduces to 

^iio_ j + 1/2 2(j-£) 
"1 - 2£+l ' "^"J " £+1/2 ' 

which is a well-known result (see e.g., [H]). Notice that, due to the conservation of jz, the product 
e''(A4)£'*(A4)* in (g]) can be replaced by S'^'^. We quote here also the simplified form of ^ in this 
case, 

5(j,AAi/2,i/2) = c},fSi,jJs'^s^ + ^vW+T)c|f (£,4;i,A4K4)(i/2,4;i,A4|i/2,s,) , (5) 

which results from ^ by applying the WE theorem. 

|A£| = 1 In this case 5(j,f 1/2, 1/2) can be non-vanishing only if £' = £ + 1 and j = i + 1/2, or if 
£' = £ — 1 and j = £ — 1/2. Equation ^ takes the simplified form, 

50\£i,£,l/2,l/2) = c||'(£i,^;|e^XA£Fl^,4)£''(A£)*XA(4)*ia%Xi3(s.) , - -2 , (6) 

which is appropriate to derive an addition theorem for spin- 1/2 spherical harmonics. As in the previous 
case, the product of standard versors can be substituted by S'^^. Similarly, ([2]) reduces to 

S{j, £,,£. 1/2, 1/2) = -A£y|^|^^(£, 4; 1, A4Ki, 4) (1/2, 4; 1, A4|l/2, s,) . (7) 

As a side remark, we notice that in this case (2(£) + l)/(2i?i + 1) (2j + l)/(2j + 1 + A£), which provides 
an alternate form for ([7]). 

2.2 s' = 1 = s 

The general result ([3]) is greatly simplified in the case of S{j,£',£, 1, 1) because, for fixed £, j 
only the three values ±1, 0. The spin matrix elements are given by 

(l,4|enA4)*5''|l,s,) --*e^A4)*(£(4r Ae(s,))'' , 
{l,s',\e '''''' {A£z)*S'''S'''\l,s,) = -£''i''^(A4)*e''U4)*e'''(s.) • 
There are three possible values for A£ — 0, 1, 2. 

|A£| = In this case ([3]) can be written as 

Sij,£,£, 1, 1) = Cg°S,,jJ,,s. - Clll{£,£',\e\A£z)L^\£,£z) £"(A4)*^ (?(<)* A e (s,))^ 
- Cllf{£, 4 1£ '^^'^^ (A4)i'=^ L''^ % 4) £ "^"^ (A4)*£ «)*£ ''^ (s.) , 

^«'° = ^ItT' ^1^' = ;^(2j + i)(0--^)0- + ^ + i) + i). C'ii2 = :^(2j + i), 

D,, = {-\y-'-^ ((j ~ £f + 1) (j + £ + 2)(j + £ + l)(j + 4 . 



— £ can take 
(8) 
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This form of ([3]) is needed in the derivation of addition theorems for vector spherical harmonics. Notice that 
in the term muhiphed by C^^^ we can : 
we can rewrite the matrix elements as 



in the term multiplied by Cj^j we can replace £^''{A£z)e'^{A£z)* by 5'^'', and in the term multiplied by C|,^^ 



The relation ([2]) in this case reduces to 



/I 
6 



/TK (11) 

+ ^VW+^VW^Wi + ^Cji'{e, 4; 2, O (1, 4; 2, A4|l, s,) , 



with C^lf" as in dH). 



I A£| = 1 In this case S{j, £' , £, 1, 1) can be non-vanishing only ii £' = £ + 1 and j = ^ + 1, ^, or if £' = £ - 1 
and j — £, £ — 1. With that discrete set of values for j — £, the spin matrix elements ([8|), and the reduced 
matrix elements in appendix [Cl from ([3]) we obtain 



4, A 1, 1) = -Q\V,(4,4|£'(A4)f 1^,4) *£^A4)*(£ (<)* A £ (s,))'' 



(12) 



As before, on the first line of ([T^ we can replace e'^(A£2)e ''(A^^)* by J'^'', and on the second line we can 
rewrite the matrix elements as 



(^',^;|£'=i'=^(A4)f'=^i'=1A4)£^^^^(A4)*£''Ks'J*e''^(s.) = (^',4l/^''i'^°K,4)e^(4)*e''(s.) . (13) 



The form of ([2]) in this case is 
S{j 



J, £,J, 1, 1) = A£^l^±^ (clllj{£,£.; 1, A4|4,^i)(l, 1, A4|l, s.) 

+\\[\^m-im£) + 2,)Clll 4 ; 2, A4 14 , 4) (1, 4 ; 2, A4 1 1, s.) 



(14) 



with the coefficients (HI' 



I A£| = 2 In this case S{j, £',£, 1, 1) can be non-vanishing only if f ' = ^ + 2 and j = £ + 1, or if f = ^ - 2 
and j = £-1 . Using dH), ([as]) and (|a3l) . from (0]) and ([2]) we obtain 

Sij, 4, ^, 1, 1) = (4,4 1£ '^'^ (A4)f '-^f \£, 4) £ '^'^ (A4)*£ '^^ (4)*?"^ (s.) 

^ ^ (2j + i)(24 + J gl.'l,(^^4;2, A4|4,0(i,4;2, A4|i,s,) , ^^^^ 

112 _ 2i 1 
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vW+i) 



The first line is the form needed for an addition theorem for vector spherical harmonics. Notice that the ma- 
trix element on that line is traceless, because £' ^ £, and symmetric, so we can replace e '^i'^^ (A^z)^ ''^''^ (A4)* 
there by Sf^^^igk^h^^ r^^^ second line is just 
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2.3 s = 3/2 = s' 

The spin matrix elements entering ([3]) when 5 = 3/2 = 5' are given by 

£'=(A4)*(3/2,5i|5^-|3/2,5,) = ^£'(A4)*xA(4)*4sxl,(5.) , 
e'^"(A4)*(3/2,5;|5* 5^13/2,5,) = -3r'^(A4)*xA(4)*xi(s.) , (16) 
e^^'=(A4)*(3/2,5;|5'5^-5'^-|3/2,5,) = -^£^-'-^(A4)*xA(4) V^^xK^.) • 
The orbital angular momentum change is < A£ < 3. 



|A£| = Substituting (HH) in © we get 

SU, £,£,3/2,3/2) = Cy°StjJs'^s. + ^C|f (^,^;|£^(A4)i1A4)£^A4)*XA(4)*4BXB(s.) 

- 3c|f (£, e',\e'^'^ {M,)V-V^ % Qe^^^^ (A4)*Xa^ (4)*Xa' (s.) (17a) 

- ^1*^*3 (A4)i'^ i^^i^^ \l, 4)£ (A4)*Xa' (4)*4^sXb (s^) . 

where the coefficients C/^J given in ([3]) in this case take the form 
|o 1 2i + 1 



3 3^r 



4 2^+ 1 



^ff3_ 1 {-xy-'-'i^ 1 

2 |j-£|(2£+l) i?,, ' 

5., - f^ + i + ^fj-^-l) fj-^ + ^l (.7-^ + 1) 



As above, we notice that the matrix elements in (jl7ap can also be written as, 

£^(A4)(^,4|i^K,4)enA4)*xA(4)VlBXs(^.) = (^,^'Ji1^,4)xA(4)*4BXl.(^.) , 
2£^''?(A4)(^,^;|i^L«|A4)£*^'(A4)*xA(4)*xM^^) = (^:^;ii^^i''^1^,4)xA(4)*xi(s.) 
6e^'«'-(A4)(A4|i''i'in^,4)r'^'(A4)*xA(4)VlBxMs.) = (^,^'.|i^*i^i'^°K,4) 

X xKOV^bxKs^) ■ 



(17b) 



(18) 



With the matrix elements written as in ([T7a|), S'(j, ^, ^, 3/2, 3/2) takes the form needed in the derivation of 
addition theorems for spin-3/2 spherical harmonics. 
Similarly, ([2]) can be rewritten as 

1, 1, 3/2, 3/2) = 4fbt'^i^bs'__s^ + 4f '(^. 4; 1. A4|^, ^;)(3/2, 5^; 1, A4|3/2, 5,) 
+ 41 '(^, 4; 2, A4K, 4) (3/2, s',; 2, A4|3/2, 5,) 

4f'(^,4;3,A4K,4)(3/2,5;;3,A4|3/2,5,) , ^^^^ 
3 ,Jii_./i~ 



33r, 33r| 33i / 15 / 3 3-, 3 3, , , 33, 

= ^_v/£(7TT)c/,| , Ki,] = ^ -yZ(ZTT) v/(2^ - 1)(2^ + 3)c/,; 



3 3 o 
.2 2'' 



2 y ^ V(^ - l)^(^ + l)(^ + 2)v/(2^ - 1)(2^ + 3)C,^ 
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Eqs. P7|) and are of course related by the WE theorem, with the reduced matrix elements in appendix 



|A£| = 1 In this case 50, f, ^, 3/2, 3/2) can be non-vanishing only if f = £ + 1 and j = £ + 3/2, ^ + 1/2, 
i - 1/2, 01 a e' = £ - 1 and j = i + 1/2, £ - 1/2, £ - 3/2. Evaluating the spin matrix elements with ^ we 
find, 

5(j- €1,^3/2,3/2) = ^a||'(£i,4r(A4)fl£,4)£''(A4)*xi(4)V%xi(s.) 

- 3cll^{£,,£'jr^^^ (A4)f ^^i'^ \£, 4) (A4)*xi^ (sz) (20a) 

3 3 3, 

2 



- :,Clj;{£^J'Je'^'^'HA£,)?^^V^V^\£,£,) e=^=^=^{M,Xx'^{s',ra%x'B'{sz) , 



3 3 A 

with the coefficients C/^l^ given by 



3 ^ _1 j + 1/2 v/(j + ffl + l)'-4v/4-(j~(€))^ 
'''' 5 2{£) + l ^ {£){{£)+!) 

(3W-J + l) (0--W)(j + l/2) + l/2)(5(j-W)2- 
^2 2^ _ 2 \ / A / V / (20b) 

m + i)V{emTT) ^ij + {i} + if-4^A-ij-{£}f ' 

Vffl((€) + 1) (3j - W + l) (3/2(j - (€))^ - l) 
= ~^2(.) + DjU + 1) ^(, + (,) + i)^_4V4-(,-(£))^ • 

This expression in terms of matrix elements of irreducible tensor operators is needed in the derivation of 
addition theorems for spin-3/2 spherical harmonics. We notice also that the matrix elements in (j20ap can 
be written without standard tensors as, 



(4,4|e'(A4)fl€,4)£''(A4)*xi(4)V%xMs.) = (4,41^1^,4) xi(4)*4BxMs.) , 
{£^,£'Je^^^^{A£,r^L^^\£,£^)e^^^^{A£X-rA'{s'Jx^^^ 

{£Me'''^'''{M.)r''l^'^L^'\^A.)e'''^''{^^^^ 

xxiH4)V:i5Xl^'(s.) . 



(21) 



Analogously, ([2|) can be written as 

5(j, 4, €, 3/2, 3/2) = 41' 1' 1^1' (3/2, 4; 1, A4 13/2, s,) 

+ 4!'(A^-; 2, A4|4,4)(3/2, s',; 2, A4|3/2, s,) (22a) 
+ 4S(^'^-; 3, A4|4,€;)(3/2, s',; 3, A4|3/2, s,) , 



with, 



3 3, 1/15 A£ , 3 3i 

= i\/f 7i|^^(2(£) - \mi) + i)m + 3) c|f; , (22b) 
ii3 1 /2r A£ 



1; = -^- ^___ V(2(£) - 2){2{e) - l){2{e) + 1)(2(£) +3)(2(£) +4)C,^|; 



33^ 

where the coefficients C^Jj are given in (j20bp . 
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|A£| = 2 In this case S{j,l' ,£,3/2) can be non-vanishing only if f = ^ + 2 and j = ^ + 3/2, £+ 1/2 or 
i' = £ — 2 and j — £ — 1/2,1 — 3/2. Taking into account this reduced set of j values, and the spin matrix 
elements (|16l) . eq. ([3]) takes the form 

Sij, £2,£, 3/2, 3/2) = -3^11/(^2, 4 {M,)r^-?^^ \£, £,)e^-^- (A4)*Xa^ (4)*Xa^ (s.) 



3 3 



2 3j--2(£) + l/2 ^||3 8 j-(^) 



2 2 ^ ± ' ' ' C' 



^'^'^ V3 V(2j-l)(2j+3) ' '^'^ V3 V(2j-l)(2j + 3) ' 
We note here that the matrix elements appearing in this equation can also be written as 

i(f , 4|f {^^f '^i'3}o 4) (4) V^^^xl? (..) , 

where the matrix element in the r.h.s. of the first equality is traceless because £' 7^ £. Similarly, from eq. ^ 
we get 

S{j, 4, £, 3/2, 3/2) = kI},^{£, 4; 2, A4|4, £',) (3/2, 4; 2, A4|3/2, s,) 

+ nlj^ie, 4; 3, A4I4, 4) (3/2, 4; 3, A4|3/2, s,) , 

/ ^112 (25) 

2^^ y {2{£)-l){2{£) + l){2{£)+3) '-'1 ' 
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ii3 1/21.;—- ((^)-i)(i?)((^) + i)((^) + 2) 33 



2V 2 V (2W-l)(2W + l)(2(^)+3) ' 
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33^ 

with the coefficients C*/^!^ of ((23|) . 

|A.£| = 3 In this case 5(j,f,£, 3/2, 3/2) can be non-vanishing only iff = ^-K3andj = £-H3/2,orf 
and j = ^ — 3/2. Eq. ([3]) then reduces to 

^(j, 4,^,3/2, 3/2) = (4,/J?"^'=^'=^(A4)f"^f"^f'=3K4) 
^^/...2^3(A4)*^;ji(4)V^^^xi^3(,^) , 



(26) 



3V (2j-l)(2j + 3) ' 



which is the form needed to formulate an addition theorem for spin-3/2 spherical harmonics with £' = £±3. 
The matrix element in (1261) can also be written 



(27) 



{£', 4|e '=i'=='=3 (A4)f '^if '==f 1£, 4)£ '^^''^''3 (A4)*Xa' (^'zT^abXe' i^z) 



= i(f ,4|f{'=^f'=^f'=3}o|^^4)^^.(,/)V^^5xlM^.) • 



Similarly, ^ reduces to 



Sij, £3, £, 3/2, 3/2) = (^, 4; 3, A4 14, 4) (3/2, 3, A4|3/2, s,) , 

333 _ 77 / 2£+l / (2ffl-l)(2(^) + l)(2(^)+3) 333 (28) 

'*^3^.- 32 ^'^y (4 + 2 V ((£) - 1)(^)((^) + 1) "^^3^^- ' 



with Cllj' from ([Ml) 
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2.4 The spin transition operator 

We define the spin transition operator T as a spin vector operator, commuting with aU orbital operators, 
and satisfying 



[S\T^] = ^e^^^T'^ , {s'\rT^\\s) = f^^M-tl (j^,^^^^^ + 5,^,^,)) . (29) 



From (123 we have, 



/ / / 1^2, \ \5s's5s'^s^ if s' > 1/2 or s > 1/2 

if.s'=.=i/2 • 

From and the operator T is seen to be the spin-space analog of the orbital operator r. In particular, 
the reduced matrix elements of tensor products of r and L given in Appendix [C] apply without changes to 
tensor products of T and S. 

2.5 s', s = 1, or 0, 1 

The relevant spin matrix element in this case is, 

(i,4|r''|o,o)^-Le''(4)* . (31) 

\A£\ can take the values and 1. 

\A£\ = S{j,£,£, 1,0) (resp. S{j,£, £,0,1)) can be non-vanishing only if j = ^ (resp. j = f). Then, © 
and ^ take the form, 

S{e,e,e,l,0) = ^Cl^,\£,£'jL%£,)e'{s',)* = -V3{i,i,;l,Al,\£,£',}{l,s',^^ , 

r-T- (32) 

^101 _ 
'^m — 

Similarly, 



£{£ + 1) 



5(AA^,0,l) = -^C'°V(^,4|L'=|^,4)e"(s.) = (^,4;l,A4|^,4>(0,0;l,A4|l,s.) , (33) 
with Cgll — C}^1. 

|A£| = 1 In this case S{j,£' ,£, s' , s) can be non-vanishing only ii j = £ (if s ~ 0) or £' (if s' — 0). Thus, 
with the spin matrix element (PT|) we get, from (jS]) and © , 

(jWl / 2/ I 1 

(34) 



/ 2^+1 
2(4 + 1 



Analogously, 

^•011 



5(4,4,^,0,1) = ^^(4, ^'Jf'=|£,4)e'=(s'J = (^,4;l,A4|4,^'J(l,s;;l,A4|0,0) , 

(35) 



I 2£' + l 
2(4 + 1 



8 



2.6 s', s = 3/2, 1/2 or 1/2, 3/2 

The spin matrix elements appearing in ^ in this case are 

(3/2,4|£"(A4)*r''|l/2,s.) = i./|£''(A4)*x^(4)*XA(s.) , 

(36) 

There are three possible values for < \A£\ < 2. 

\A£\ = lis' = 3/2, s = 1/2, andf = £, then S'(j, 3/2, 1/2) can be non-vanishing only if j = ^±1/2. 
In the case s' = 1/2, s = 3/2, it must be j — £' zL 1/2. Substituting (|36|) in ([3]) we obtain 



^0\^,^,3/2,l/2) = iy|c|f (£,fj£^A4)2.'=K,4)e''(A4)*x5(4)*XA(s.) 



(37a) 



with 



^f^i_ 2 2j-^+l/2 ^H2_ J-^ 1 

- 2£+lV2^-J + l/2 ' """^ ^^£+l^(2j-£+l/2)(2£-j + l/2) ' ^ ^ 

The matrix elements in p7ap can also be written without standard tensors as, 

{£J'Je\Ae,)L%£,)e'^{A£,rXA{'^',rXA{s.) = {£J'jL%£,)xXi.<rXA{s.) , 
{£, 4 1?'=^'=^ (A4)i'=^i'=^ \£, 4)£ "^"^ (A4)*X^^ (4)*<^abXb(s.) = 

= i(^,4|Lt^-^L'=^>"|^,4)XAn4)*4^5XB(..) . 

(38) 

Similarly, ^ takes the simpler form, 

S{j, £, £, 3/2, 1/2) = - ^ yZ(ZTT)c|| ' (A 4 ; 1, A4 £',) (3/2, 4; 1, A4 1 1/2, s,) 

(39) 

-^VW+1)V(2^- l)(2^ + 3)c|f (^,4; 2, A4|^, ^;)(3/2, 4; 2, A4|l/2, s,) , 

with the coefficients given in p7b| . 

The case s' — 1/2, s = 3/2 is completely analogous. From wc find 

5(j- 1/2,3/2) = i^/|c|f(A4|e^A4)Ll^,4)£''(A4)*XA(4)*xi(s.) 

2V2^ (40) 

1 1 3 

IO^--2,. Di\^ktk2(Ar, \TkiTk2\f \-S-hih2/An .( J \*h2 ^hi 



with C^^J = C'; (see (l37b|), and from © 

5(j, £, £, 1/2, 3/2) = \ ./|y^(ITl)c|f (€, 4; 1, A4K, 4) (1/2, 4; 1, A4 13/2, s.) 

^ (41) 
1 /5 



+ g y 2 ^^^^+^^"(2^ - 1)(2£ + 3)C/,f (£,4; 2, A4K, 4) (1/2, 4; 2, A4|3/2, s,) . 

The expressions ([37| and (|40|) will be used below to derive addition theorems involving one spin-3/2 and 
one spin- 1/2 spherical harmonic. 
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I A£| = 1 We consider the case s' — 3/2, s = 1/2 first. In this case S{j, i', £, 3/2, 1/2) can be non-vanishing 
only if f = ^ + 1 and j = £ ± 1/2, or if £' = ^ - 1 and j ^ £±1/2. With the spin matrix elements in 
from ^ we obtain 

SiJJ^,i,3/2^/2)^l^^lcl|;{e^e',\e''{Ae.)r%e.)e^^ 

Y (42a) 

with the coefficients 



; (42b) 

^ _ j_ /I + AsAeU - {£)) 1 + 2A,sA^(j - {£)) 



where in this case As = 1. On the first line of (|42ap we can replace e''{A£z)£'^{A£z)* with 6''^, and on the 
second line, 

(/f,|£'=^'=^(A4)f'-^L'=^K,4)e''-^''-^(A4)*xl^(4)V!^^BXB(s.) = 

= l{£'£'Ari'^L>^^^°\£,£.)x!l'{s',ra%XBis,^) . ^^^^ 

Similarly, from ^ we get 



SU, £,,£, 3/2, 1/2) = -l^H^^ ^TT^^i'^^'^- 1' A4|4,4)(3/2, 1, A4|l/2, s,) 



(44) 



30^^./2(£) + l 



32 



A^^(2(£)-l)(2(^)+3)C,^|;(^,4;2,A4Ki,4)(3/2,s;;2,A4|l/2,s,) 



with the coefficients ()42bp . 

In the case s' = 1/2, s = 3/2, 5(j, f , ^, 1/2, 3/2) can be non-vanishing only ii £' = £ + I and j = £ + 3/2, 
^ + 1/2, or if = £ — 1 and j = £ — 1/2, € — 3/2. The treatment of this case is completely analogous to the 
previous one. In terms of matrix elements, from ([3]) we obtain 



^(j;£i, £,1/2,3/2) = ly|c|J(44|e'XA4)f'=|^/.)?^A4)*XA(4)*XA(s.) 

1 1^ 1 3 



(45) 



+ -^\l7fU;{h£;A^''''''{A£;)r^-L^^%£^)e>^-'^^ , 
with the coefficients C'/J^ given by (I42bp with As = — 1. From ([2|) we get, analogously, 



^(j, 4,^', 1/2, 3/2) - ^A£^l^±^Cll^ {£,£,; 1, A4|4, 4)(l/2, s^; 1, A4|3/2, s,) 

+ ^A£^^^±iv/(2(^) - 1)(2(4 + 3)c|f;(£,4; 2, A4 14, 4) (1/2, 4; 2, A4|3/2, s. 



(46) 



The equalities (jH]) and (|^5|) will be used for the derivation of an addition theorem involving one spin-3/2 
and one spin-1/2 spherical harmonic, with orbital angular momentum differing by one unit. 



10 



|A£| — 2 As in the previous paragraph, we consider the case s' = 3/2, s — 1/2 first. In this case 
S{j, e'J, 3/2, 1/2) can be non-vanishing only if £' = £ + 2 and j = £ + 1/2, or if f = ^ - 2 and j = £ - 1 /2. 
Using the spin matrix elements ([55]) and the range of j, ([3]) reduces to 

1 /3^|^2,^ e'\p-klk2/Ae \C:-kic-k2\n n \^hih2 { A n \* c'^i t J \* 



S(j,f2, £,3/2^/2) = JW2CfJ;(44|e"'"'(A^.)F'■f"'|t^.)£»■»■(A4)^y■(»l)Vy^XB(..) . (47a) 



c|f;.^A.Af^^ji|±l^, (47b) 

where As = 1. The matrix elements in (|47ap can be written without standard tensors as, 

Similarly, in this case ^ takes the simplified form 



S(, £. £ 3/2 1/2) - -y^VWTl. mW+l) c§|2 

b{j,£2,£,3/2,l/2)- g (2(^)-l)(2(^) + l)(2(^)+3)^^^^^- (49) 

X (^,4;2,A4K2,4)(3/2,s^;2,A4|l/2,,s,) . 

In the case s' = 1/2, s = 3/2, S{j,£',£, 1/2,3/2) can be non-vanishing only if f = £ + 2 and j = £' - 1/2, 
or if £' — £ — 2 and j = ^' + 1/2. In this case reduces to 

with C/^l^ given by the same expression (j47bp as C^^jj , with As = —1. The equality ([2]) takes the simplified 
form 



S{j,£„£,l/2,3/2) = - J-C,^|,^(44|£^-^^-^(A4)f'=^f'=1^,4)e''^'^nA4)*XA(4)V!i^5XBHs.) , (50) 



139 

»2 2^ 



5(j, 4, A 1/2, 3/2) = — ^/2£Tly (2(4 - 1)(2(4 + 1)(2(4 + 3) ^^^^^^ (51) 
X (£,4;2,A4|4,4)(l/2,si;2,A4|3/2,s,) , 

3 Angular momentum projection operators 

We define the angular momentum projection operators as 

e 

Vt^ \£,m){£,m\ , < £ < 00 . (52) 

It is immediate that Vi are a sequence of orthogonal projectors, which commute with L and resolve the 
identity, 

00 

VtVt = VeSe'e , ^ = 1 , [L,Ve] = . (53) 



The matrix elements of Ve in configuration space follow from their definition (|52p and the addition theorem 
for spherical harmonics 

V,\?) - J2 yUrmm) , (rW) = ^^P,{?.?') . (54) 

rn— — i 

From ([52|) -([54 l) the matrix elements of Ve with other orbital operators can be computed. Those matrix 
elements, together with the results on products of CG coefficients of the previous section, are the basic 
building blocks needed in our derivation of addition theorems for spin spherical harmonics. As will be 
discussed in II, however, the results for Ve matrix elements constitute by themselves addition theorems for 
spherical harmonics. 
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3.1 Matrix elements with tensor powers of L 

The matrix elements of with tensor powers of L will be used in the derivation of addition theorems for 
spin spherical harmonics. For the irreducible components, from which the full tensor matrix elements can 
be reconstructed, we have 



{?'\L^^^ . . . L'^^^^VA?) = ^^^iP (f A V)^''" . . . (f A Vf'""' PAx) 

An 

= i (f A V)^''" (f . . . L''''-1>»-P^|f ) (55) 
(r A V)^'''' (f iL^'^i . . . L'''-^^°^°Pe\r) . 



ip-iy. 

Thus, we can obtain these irreducible matrix elements recursively, with the result 

47r 



with ' the fc"^ derivative of P^ and, 

x = r-r', v = rhr', C^,, = (2g - 1)!! (^^^J , Cp,o = 1 , ^^^^^ 

2-ii...ir.i;ji...J„2 = f *i . . .f^-if ^ _ f ina 

For notational simplicity in what follows we adopt the following useful conventions. In - yi - J"2 one 
or both sets of indices may be empty, 

Z = f , _ ^ ^ . . ^ Z ' = 1 . (56c) 

In the term g = in (|56aP the index sets are empty, ^'ii - 'io:/ii...'io ^ z ' . Similarly, in the term q ~ [p/S], 
when p is even, we set v'^^g+i _ ^ _ y'*? = yh-p+i _ ^ _ ^ftp ^ 

The cases n = 1, 2, 3 will be needed for the derivation of addition theorems for spherical harmonics of 
spin 1/2, 1 and 3/2 so we quote them here explicitly, 

{r'\LVi\r) = i^-^vP^i? ■ f ) , (57a) 
An 



(r'\Li'L^^°Ve\?) = — ( v^'v^^'> P;'(r ■ f ') + ?^'?'^^°P;(r ■ r')) , (57b) 

47r v / 

{?'\L^'L^L''^«Ve\r) = -i'^^^ (vl'v^v''^'>P^"ir ■ r') + 3f ''fVJw'=>«P;'(f • r')) . (57c) 

47r V / 



3.2 Matrix elements with tensor powers of r 

The matrix element {p'\Vi'r^^ . . .r'''"Vi\p) is an irreducible tensor ii £' — £ ±n, otherwise reducible. The 
irreducible matrix elements will be used in the derivation of addition theorems for spin spherical harmonics. 
They satisfy the recursion relation. 



1 



2£„ + 1 



{£n+l - £n)^y - -{in+1 + in + l)p" 

with 



{p\Vi„?''...?'"Vi\p), (58) 



{p'\ri,?Vi\p) = ^ {p'Pi (x) - pP^ix)) . (59) 
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In these equations we used the notation = izLk, see appendix El Eqs. (1551) . ([5^ are established in 
appendix |D] The recursion relation ((55)) can be solved with the initial condition ((5^ to give 

fci+fe2=" (60) 

47r nfe;o(24 + l) 47r^ n J ^ + i + n - ^y.l' 

with ^*i- -'"iyi -J"2 defined in (j56bp . (j56cl) . As before, we adopt the convention that in the term with ki — n, 
A:2 = in (|60|) the first index set in Z is empty, ^*i - «o;*i - «" = z ;'i - *'> and, in the term with ki = 0, /c2 = n, 
the second index set is empty ^*i - «";*i>+i - *ii = ^ii...i„; ^ 

The particular cases n = 1, 2, 3 will be needed for the derivation of addition theorems for spherical 
harmonics of spin 1/2, 1 and 3/2, so we give them here explicitly, 

{p'W.r-Pelp) = ^ i-p'Pii^) +p"PLi^)) , (61a) 
{p'\V,,r'r^-V,\p) ^ Qp{''^p'^>°P;'(a;) -p{^^p'*^>"P;:(x) + ^p'i^^p'-^}op^'^{x)^ , (61b) 

{p'\re.r'r'^r^"Pe\p) = ^ ^^^^ ^ ^ (^-^p^^^p^^p^^^» {x) + ip^'^p^^p "^>"P-(:r) 



(61c) 



3.3 Mixed matrix elements 



We consider now matrix elements of Ve with tensor products of both L and r . The matrix elements of the 
irreducible tensor operator r^^^ . . . r '^^L'^^ . . . L'^p}o are of the form 

(^'\ViS^'' . . .P"L'*i . . . L'''^"Vi\p) ^i\p A Vp)^''* . . . (p A Vp)'''{p'\ViJ'' . . .P">°P^|p) . (62) 

From this equation and (IM)) . for t = 1 we obtain 



{p'\ViJ^'' ...r'"L''^''Vi\p) ^iA^ J2 (-l)'=^('^~)z{*i-'^-2'"=2+i-'"z;''>«p/"+^^(a;) . 

fcl,fc2=0 ^ 

ki+k2—n 

For s > 1, equation (|62|) leads to the recursion relation 
{p'\Ve„r^'' L^' ... L^'^^Velp) = 

= , , ' ,„ (P A Vp)f"(p'|P,„f . . .f ^"i'^i . . . i'''-i>«>"P,|p) 
(n + i — 1)! 

which can be solved with the initial condition ((63|) to yield 



(63) 



(64) 



(65) 



ki,k2=0 ^ ' 9=0 

fcl+fc2=" 



with [ ] in the inner summation denoting integer part, A„ as defined in (j60p . and Cs.gi 2;, iT, and ' yi - J"2 
as defined in (j56bl) . with p instead of r . As before, for notational simplicity, we have not separated from 
the sum the terms with (k\.,k-2) = (n, 0), (0,n), those with 5 = and, for even s, those with q ~ [s/2]. In 
those cases we apply the same conventions as explained after eqs. ([55)1 and ([50)1 . 
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Particular cases of importance for the derivation of addition theorems for spherical harmonics of spin 
1/2, 1 and 3/2, are 

{p'\V^,?^^L^^V,\p) = {p'^^v'^^P',[{x)-p^^v^^P'/{x)) , (66a) 

l$'\Vi,9^'?^L''^«Vi\p) = ^777^ (p't^p P;;'(a;) - 2p'{''p%'=>''P;;'(x) +p^'p^v''^«P'i/'{x)] , (66b) 

47r zt]^ 1 V / 

(p'|Pf,f{'L^'iL^'^>°-p<;|p) = --^^ (p'^*z;^'it'^'^>°F;^'(x) -p{*wJ'it'^'^>°P;''(a;) 

-p'{^p^pi^}^'P'/{x)^ . (66c) 

4 Final remarks 

We have presented in the foregoing sections preliminary results needed for the systematic derivation of 
addition theorems for spin spherical harmonics in II. In sect. [2] we obtained the factorization of orbital and 
spin degrees of freedom in products of CG coefficients of the form ([T]) in general form, and discussed the 
particular cases with < s',s < 3/2. In those cases the coefhcients Cf,l^, given in full generality in ([3]), 
were reduced to much smaller forms, and the tensor and spinor structures of each term in the expansion 
given explicitly. 

In sect. [3] the matrix elements of the angular- momentum projector operator with the irreducible compo- 
nents of arbitrary tensor products of r and L are given in general form, for all values of their parameters. 
Those matrix elements will be used in II to obtain general expressions for bilocal spherical harmonics, and 
to derive addition theorems for spin spherical harmonics. Their applicability is even wider, however, when 
they are appropriately combined to obtain matrix elements of reducible tensor operators. Some examples of 
those applications will also be considered in II. 

As a side remark we point out that an unexpected byproduct of the results of sect. [5] is an improvement 
in computational efficiency, at least in the specific case of infinite-precision computation |13j in which the 
results are given as a rational number times the square root of a rational number. Consider, for example, 
the computation of both sides of p^ . at fixed (. and j, for all possible values of < i'^Az < ^ and 
—3/2 < s^,Sz < 3/2. If the l.h.s. of (IT9)) is computed as written on the second line of ([T|), we find that 
the ratio of CPU times Tr.h.s./'T'i.h.s. begins at ~ 2 at ^ = 1, monotonically decreasing with i to reach 1 at 
i 20, ^ 0.5 a.t £ ^ 100, and ^ 0.25 a.t £ ^ 200. We remark that those numbers are subject to statistical 
fiuctuations. 
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A Notation and conventions 

Throughout the paper we set ?i = 1. We denote tensor product states of orbital and spin angular momentum 
by \(-Az]s,Sz) = \^,^z)®\s,Sz), and states coupled to total angular momentum j, by \(-,s,j,jz)- We follow 
the notation of [7] for CG coefHcients, 

{iAz;s,sz\j,3z) = {e,ez;s,sz\e,s,j,jz) ■ 

(Different notations are used, e.g., in [UinilHlIIl]-) We adopt the usual convention that {£,£z;s,Sz\j,jz) — 
ii iz + Sz jz, or j < 1^ — s|, or j > £ + s. For Wigner 6j-symbols [l]|4]-[6] we use the definition and notation 
of [HE] . The most important property of 6j-synibols for our purposes is the relation [I] 

32+ ja 

{ji,mi;j2,m2\j',mi +m2>(j',TOi + j3,m3\j,m) = ^ V(2j' + l)(2j" + l)(-l)Ji+J2+j3+j 

J" = |j2-J3| (A.l) 

X I •'^ ^„ I (ji,mi; j", 7712 + m3|j,m)(j2, "T,2;j3,"T,3|j", 7712 +m3). 

Jo J J J 

This relation is satisfied by the implementation of CG coefficients and 6j-symbols in the software system |13) . 
Our notation for the reduced matrix elements entering the WE theorem is explained in appendix IB. II For 
spherical harmonics Yim and Legendre polynomials Pi we use the standard definitions [l]|6l[8] . Furthermore, 
we adopt the convention, usual in numerical computations [13j, Ygmir) = if \m\ > i. When states with 
different angular momenta or spins are considered, we use the notations, 

= ^ , M^e-i, Mz^^z-iz. (A.2) 

and similarly (s). As, As^. We also find useful the notation — £ ± k for orbital angular- momentum 
quantum numbers, with the convention that when several , ■ ■ . , (kp appear in the same equation, then 
either the upper or the lower sign is chosen for all of them simultaneously. Thus, £n+i — (n = £p+i — (p = ±1 
for any n and p. 

We denote tensor indices by lower-case latin superindices, ■■■»'> ^ and spinor indices by upper-case latin 
subindices, XA- Since we consider only tensors and spinors of the su(2) algebra (as opposed to the sl(2) 
algebra), there is no need to raise spinor indices or lower tensor ones. We denote normal spinors and tensors 
by a caret, e.g., r — r/r. Our choice of orthonormal bases for spinors and tensors is explained in appendix IB] 

Given a numeric or operator tensor we denote its associated totally symmetrized and antisym- 

metrized tensors as, 

= ^ A*"! , ^['1 ■ ^sgn((T)A*-i , (A.3) 

a a 

where the sums extend over all permutations cr of ii, ... , i„. The same notations apply to tensor products, 
e.g., r^'^L^^ — r^U -\- U. We denote the traceless part of A'"^ - *" by A(*i- -'")o. The cases of importance 
in this paper are rank-2 and 3 tensors, 

A(ij)o _ Aij _ }_AkkHj 

3 

2 

^(pgr)o _ j^pqr _ _IJl^pjj^qr _|_ ^iqi gpr _|_ J^iir ^pq\ (A. 4) 

5 

_|_ }_^J^qjj gpr j^rjjgpq ^ y^ipi gqr _^ j^iri gpq ^ j^Hp gqr _|_ j^iiqgpr-^ 
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The same notation applies to tensor products, such as r^^r^^° = r^r^ ~ 1/3(5*^ . We denote by the 
traceless part of y4'f*i - '"h Thus, the irreducible component of A*^ - '" is - . A more systematic 

approach to irreducible tensors is given in the following section. 

B Standard tensor and spinor bases 

In this appendix we introduce the standard bases for irreducible spinors and tensors of the su(2) algebra used 
as spin wavefunctions throughout the paper. We use spin wavefunctions of the Rarita-Schwinger type 14 , 
which carry only vector indices varying from 1 to 3 and spinor indices from 1 to 2, and whose components 
are linearly related by symmetry and tracelessness and, for spinors, by transversality conditions. In fact, 
the spinors discussed in sect. IB.2] are the non-relativistic version of Rarita-Schwinger spinors. The standard 
bases are defined so that they are orthonormal, have definite complex-conjugation properties and satisfy the 
Condon-Shortley phase conventions [T1I3|31[5]. 

Our starting point is the usual basis for Pauli spinors. 



Xa(1/2)= (^J) , xa(-1/2)= (^J) 



(B.l) 



The triplet component of basis-spinors tensor products yields a basis for the representation space of the 
adjoint representation. 



XAB{m)= ^ {^,s^;^,s^\l,m)x{s'^)AeBcx{sz)c , m= 1,2,3, (B.2) 
whose associated unit vectors are. 



These are the standard basis vectors, satisfying, 

1 

s'is^r = i-iy^s'^i-s,) , sHs.ye\s',) = <5,,,, , ^ e\s.)e'^' {s,)* = j'^'^' . (B.4) 

Furthermore, from the action of the spin operator on basis spinors we get 

e'M^BXsGs,) = ^J^Xa((-1)™s,) , m = -1,0,1,5, =±1/2, (B.5) 
Vl + 

a relation that will be needed in the discussion of spinors below. 
B.l Standard irreducible tensor bases 

Starting with the standard vectors (jB.3P we define the standard basis of irreducible rank-n tensors recursively 
as, 

1 ri-l 

e'i---*"(TO) = ^ ^ (n — 1, to'; 1, s^jn, m)e'^'"*""^(m')e*"(s2) , ~n < m < n . (B.6) 

Sz— — 1 m' — — n+1 

This recursion can be solved to yield, 

1 

e*--"(m)= fn{si,...,SnW'{si)...e'-{sr,) , (B.7a) 



Si ,...,S^i — — 1 

si+...+Sji— m 



with 



...,.„) = n(i, ,s, ; , - 1, g s.\j, g = (^(^ Y{ia^ + s,)Mi-s,y. ) ■ (^-^^^ 



16 



The second equality in (IB.7b[) follows from the first one by using the explicit expression for CG coefficients 
coupling angular momenta differing in one unit [HIH]- From (IB. 71) it is clear that e'^ - *" (m) is totally 
symmetric. Explicit evaluation shows that e*i'^(m), —2 < m < 2, arc traceless and therefore it follows from 
(IB.6P by induction that e 'i- -'"(m) are totally traceless for all n > 2, —n < m < n. Thus, (jB.6p defines a 
set of 2ri + 1 irreducible tensors of rank n. From ()B.4|) and (|B.6p we obtain by induction the orthonormality 
and complex-conjugation relations 



(B.8) 



The completeness relation for the basis of standard tensors defines the orthogonal projector x*^ --*"'^^ --'" in 
the space of rank-n tensors onto the subspace of completely symmetric and traceless tensors 



X 



ll...ln;ji...]ri 



" 1 



(B.9a) 



Here, the prime on the r.h.s. denotes symmetrization over the indices zi . . . i„ and ji ■ ■ -jn separately, for 
example. 



j^ni2;jij2 _ 1 I (5nji^i2j2 _|_ jiU2ji2ji _ ~^n»2^jij2 



- E 

6 ^ 



cr{l,2,3} 



5 10 10 



(B.9b) 



jji2i3 



5 10 10 

5 10 10 



where in the last equality cr{l,2,3} is the set of permutations of {1,2,3}. Using the recoupling identity 
(jA.ip and (jB.6[) we can derive the important relation. 



n — q 



*"(s^)= ^ ^ (n-q,m';g,s^|n,m)£'i-*— '(mOe'—'+'-'-Cs^) , 1 < <? < n - 1 , (B.IO) 



Sz— — q m' — — n+q 



showing that the maximal-rank coupling of two standard tensors is again a standard basis tensor. From 
(|B.7|) we obtain the equality 

ri-*"(±n) =rH±l)---e'"(±l) , (B.ll) 
which is useful to evaluate reduced matrix elements. We adopt the convention 



£^-'"(s,) = if \.s,\>n. 
The matrix S'(„) of the spin operator 5 in the basis of spin-n states with n integer is given by 



(B.12) 



_ j>ii...i„ Qi ...i„:fei...fc„^fci...fe„ ^ 



j \ii...ir,;ki...k„ 



(B.13) 



9=1 



p=l 



If is a rank-n real tensor, or a self-adjoint tensor operator, its irreducible component is 

n n 
n\ ^ — ^ ^ — ^ 



Sz — — n 



Sz——n 



Furthermore, for self-adjoint tensor operators the WE theorem ^ holds, 



(B.15) 
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For a normal vector r, let Xim(^^) = e ' (TO)r *i . . .r*". It follows from (jB.6l) by induction that 

L^ynmir) = n{n + l)3^„™(f) and L^3^„™(f) = m3^„„(f) . (B.16) 

We omit the proofs for brevity. We must then have ynm{r) oc Ynm{r). Evaluation of both sides of the 
proportionality relation at r = z yields 



r,„(f) = -L^^I^I±^e^^-^^(m)?^^...r^^ . (B.17) 



20F 

From this relation and the addition theorem for spherical harmonics we get 



J2 e''-"'{m)e''-^'{r 
^=-e (B.18) 



(2^-1)!! 

The relation (jB.17[) can be inverted, for each £, by extending the basis (|B.6p with appropriate reducible 
tensors to an orthogonal basis of the entire space of rank-n tensors. 

B.2 Standard spinor bases 

We define the standard spin-(n +1/2) spinors with n > 1 as 

1/2 n 

xT-""{sz)^ J2 51 <^'™;l/2,4l«+l/2,s,)£*-'"(m)xA(4) , -n-l/2<.s, <ri + l/2 . 

s'^=-l/2 m=-Ti 

(B.19) 

Using (jB.6p and the recoupling identity (lA.ll) . from (jB.19l) we find the equivalent expressions, 

g+l/2 n-q 

J2 E ("-'7,"^;9 + l/2,4l" + l/2,s,)e^^-^"-''(m)xl"-''+"-'"(s;) , (B.20) 

s'^= — q— 1/2 ■m=-n+q 

with 1 < q < n — 1. From definition (jB.19[) we see that xA^ completely symmetric and traceless in its 

tensor indices. Thus, (|B.20I) remains valid after an arbitrary permutation of tensor indices in its r.h.s. From 
(|RT9| with n = 1 and (|R5|) we find, for the spin-3/2 spinor, a^gxhisz) = 0, -3/2 < < 3/2. Thus, from 
(|B.20p with q = 1 we immediately obtain, 

(^a'sXb (s.) =0 , l<k<n, -n-l/2<s, <n + l/2. (B.21) 

From (IB.19P and (jB.21[) . the number of independent components in the spinor basis is 27z + 2 as it should. 
The orthornormality and complex-conjugation relations read 

xT-'"i^zrxr-'"i<) = Ss.s'^ : <Bxh'-'"{^z) - {-i)'^'^''xr-'"i-^zr ■ (B.22) 

The completeness relation for spin-s spinors, with s ^ n + 1/2, n > 1, defines the orthogonal projector 
^*i^.»,.ui...j„ ^^^^ ^j^g subspace of spin-(n + 1/2) spinors of c(2"+i)x2^ 

(B.23) 

r I 1 1 



"/T- ^ ckihi i;k„h„ c , ^ ■ kirhi r srk-zh-z xk„hr. 
2^1^ ■•■^ ^^^^'i^l'' ^^^^ ---^ 



J^^l ■■■in',ki...kn 



■^For notational simplicity we identify the space of complex rank-n completely symmetric and traceless tensors with c'^"+^' , 
which has the same dimension. 
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where - - is the projector defined in (jB.9p . The orthogonal complement of the spin-(n + 1/2) 
subspace in ([;;(2n+i)x2 jg ]ier(X), which can be parameterized as ct^b^b with - spin-(n— 1/2) 

spinor. This is the spin-(n — 1/2) subspace of c'^"+^)^^. 

The spin operator for s = (n + 1/2) spinors is given by the operator tensor product of the spin-n and 
spin-1/2 operators, its matrix being 

iSl+./^^yiB'"-'''---'" - i^^^^^^^ , (B.24) 

with S'^^j defined in (|B.13p . From (jB.21[) and (IB.24P we obtain the following useful identities, 

ie'''''^a'XBXB'-'-{sz) = xr-'"{s.) , (B.25) 
4bXb"-'"(^.) - ^abXb^-^"(^.) + ^e'^^'^x'^■■■^"i^.) = , (B.26) 

= (n + l/2)x^-''"-^(4)*^eP'=«Xs"■■'"-^'(s.) , 
which remain valid after permutation of ii . . . in, etc. 



C Reduced matrix elements 

In this appendix we gather expressions for matrix elements of certain irreducible tensor operators needed in 
the foregoing. The notation used is the same as in (jB.15[) . We begin from the known results [8] 



{fWe'J'Wj) - ^/W+T)6,,, , (C.l) 

{i'mm = -^^[^^^§^^{e,0;L,0\i',0) . (C.2) 

In (jC.ip J* can of course be any angular momentum operator: orbital L*, spin 5* or total J*. For tensor 
products of angular momentum operators we have, using (|B.11I) and (jC.ip . 



tflir.......... = jL_ J_^ ^^2l±^s„, . ,c.3, 

For convenience, we quote here the explicit form of the CG coefficient appearing in (jC.2p . If t is odd the 
CG coefficient vanishes, and if t > is even 



(^,0;n + i,0|€„,0) = (-l)*/2 



_ „ / (2n + t-l)!!(t^ r—— l ijin+i + n + tyiy. 
^ ' (n + </2)!(i/2)! ^''''^'^ {{in+£-n-t)/2)l 



X \ / t:^ J 7TTT I t> even. 

Y {en + e + n + t + iy.\ ' 

Next, we consider tensor products of the position versor r = f/r. From (|B.17p and (jC.2p 



(C.4) 



In the case |A£| = n, which is of interest to us, a more explicit expression is given by. 



■■ 2«/2 y2((£) + §) + i\/r((^) + f + i) 



/ r(W + t + i) 
r {{£) - § + 1) 



(C.6) 
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The cases n = 1, 2, 3 are of particular importance to us. Noticing that (—1) 2 — ( JSTI j ' 



Equations (IC.6P and (jC.7p remain valid under the replacements ^ — > s, — > T', that give the reduced 
matrix elements (s ± n||e*i - *"r*i . . . T'" ||s) of tensor products of the spin transition operator for integer or 
half- integer s. 

For tensor products of r and L, using (|C.3p . (|C.6p and (IB. lip we obtain, 



V2£+T / (2(£) +n+l)!r(ffl - |A^|/2 + l/2) 
2{e) + \Ae\ + iy (2(^)-n)! r((^) + |A£|/2 + l/2) 



(C.8) 



This equation reduces to (jC.Sp for Ai? = 0, and to (jC.6p for |Ai?| = rt. The particular cases n = 2, 3 of (jC.Sp 
are important in the foregoing, 

Ae 1 



{e,\\e^^?^Um = — ^___ v/(2ffl - l)(2ffl + l)(2ffl + 3) , (C.9a) 
A£ 1 



{h I r^'r^L^L'^ 1 = VW+T V(2(^> - 2){2{£) 1){2{£) + 1)(2(^> + 3)(2(^) + 4) , (C.9b) 

^^^11'^ V3'^'' + '\/(2W-l)(2W + l)(2W+3) • ^^-^'^ 

The following mixed reduced matrix elements are also useful when dealing with tensor decompositions in 
irreducible components. 



mri?A Ly\\i) - ((^1 + ^ + 1) - 2A£) , (C.lOa) 



A X>||.) = -1 - .)(2£, + 1) - 3) ^^^^A|±1L_ , (c.lOc) 

With the replacements i ^ s, Ai As, f' T\ U S\ ([OSl and ([CJ|) give the reduced matrix 
elements (s ± n\ |eH - >*"T*i . . . T^i^-i S'^ia-i+i ... 5'" | |s) for integer or half-integer s, and similarly for (jClOp . 

D The matrix elements (p'|7^£^r . . . r *"7^^|p) 

In this appendix we derive ()58p and (|59p. First, we consider the matrix element. 



J d\P,Ap'-q)p'-qP,{p-q) (D.l) 

^^($'\v,,^,v,\p) + ^±-^{p'\v,,^,r,\p) , 
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where the integrals extend over the unit sphere, and in the last equality we used the relation 

Applying the projector property of Vt in the last line of (jP.ip we obtain, 

lS'\V,,p' ■ rVe\p) = + e + ■ (D.2) 

Similarly, using 

Wp'p' ■q=p'A{qAp') , 

we get, 

{p'lVe.p' A{r Ap')Ve\p)^W,, (^^L_ip,^ . p')^ = _^L_ip ' a (p A • p') . (D.3) 
Summing (|D.2p and (jP.Sp we obtain (|59l) with the help of the relation 

xPiix) - (n - fc)P,,(x) + - ^„)(^„+i + (n + l)Pi, = Pi^, , < A: < n . (D.4) 
The derivation of (|58|) runs along the same lines, 

= j d^q{p'\Vt^^,p''-^'p' ■?ViM{q\Vt?'' ...r'-ri\p) (D.5) 
+ J d^q{p'\Ve„^,i?'-+' ~p"-+'p' ■r)V,Jq){q\Vi?'' ...r'-Vi\p) . 
Applying now (|D.2p to the first integral and (jP.SP to the second one in the last equality, we obtain (|58|). 
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